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Abstract 

We derive a closed expression for the SU(2) Born-Infeld action with the sym- 
metrized trace for static spherically symmetric purely magnetic configurations. The 
lagrangian is obtained in terms of elementary functions. Using it, we investigate 
glueball solutions to the flat space NBI theory and their self-gravitating counter- 
parts. Such solutions, found previously in the NBI model with the 'square root - 
ordinary trace' lagrangian, are shown to persist in the theory with the symmetrized 
trace lagrangian as well. Although the symmetrized trace NBI equations differ sub- 
stantially from those of the theory with the ordinary trace, a qualitative picture of 
glueballs remains essentially the same. Gravity further reduces the difference be- 
tween solutions in these two models, and, for sufficiently large values of the effective 
gravitational coupling, solutions tends to the same limiting form. The black holes 
in the NBI theory with the symmetrized trace are also discussed. 
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1 Introduction 



Recent development in the superstring theory suggests that the low-energy dynamics of 
N coincident Z)-branes is described by the SU(N) Yang-Mills theory governed by the 
Born-Infeld type action JT[|. A precise definition of such non-Abelian Born-Infeld (NBI) 
action was the subject of a vivid discussion during past few years 0, |J ^ (7[, for an 
early discussion see ||. An ambiguity is encoded in specifying the trace operation over 
the gauge group generators. Formally a number of possibilities can be envisaged. Starting 
with the determinant form of the U(l) Dirac-Born-Infeld action 
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S = — J (1 - V" det(<k„ + d*x (1) 

one can use the usual trace, the symmetrized or antisymmetrized 0] trace, a calculation 
of the determinant both with respect to Lorentz and the gauge matrix indices ||. Alter- 
natively one can start with the 'square root' form, which is most easily derived from ([!]) 
using the identities 

det(57u, + F^) = det(g^ u - F^ u ) = det(g^ + iF^) = (2) 
det(^ - iF^) = [det((^ - F%){g lw + F^))} V * , (3) 

where F 2 V = F^ a F a u (similarly for F^), and 

F F a — p p a = -n F rF o13 

F m F a v = —g^F^F^, (4) 



leading to the equality 




det(V„ + F^) = J-det(g) \1 + -F* - — (FF) 2 , (5) 



with F 2 = F^F^ , FF = F ilv F ilv . For a non-Abelian gauge group the relations (ffl) 
are no longer valid, so there is no direct relationship between the 'determinant' and the 
'square root' form of the lagrangian. Therefore the latter can be chosen as an independent 
starting point for a non-Abelian generalization. There is, however, a particular trace 
operation - symmetrized trace - under which generators commute with each other and 
therefore both forms of the lagrangian remain equivalent. This definition is favored by 
the adiabaticity argument, as was clarified by Tseytlin ||. Restricting the validity of the 
effective action by the no-derivative approximation, in the non-Abelian case one has to 
drop the commutators of the matrix-valued F^ v since they can be reexpressed through 
the derivatives of F^. This corresponds to the following definition of the action 

1 



S = — Str / \ 1 - V- det(^ + FfJ) \ d*x (6) 
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where symmetrization applies to the field strength (not to potentials 0). This action 
reproduces an exact string theory result for non-Abelian gauge fields up to a' 2 order. 
Although there is no reason to believe that this will be true in higher orders in a', the Str 
action is an interesting model providing a minimal generalization of the Abelian action 

[0- 

Some general technique was developed to deal with the symmetrized products of 
gauge generators in a symbolic form, but for many purposes it is more desirable to have 
the lagrangian explicitly. In the general case an evaluation of the action can be performed 
through an expansion of the square root in powers of F^, then the symmetrized trace over 
generators can be computed explicitly. The next step, resummation of the series into a 
closed expression, is very problematic. However this can be done if we restrict the field by 
some particular configurations (the ansatz has to be consistent with the full equations of 
motion). One example of this kind is the computation of an explicit action for DO-branes 
in three dimensions ||. Here we present such a calculation in the four-dimensional SU(2) 
theory restricting the field configurations by requirements of the spherical symmetry and 
staticity. Such configurations are encountered in the study of magnetic monopoles and 
other solitons in the NBI theory (we use a term 'soliton' in a wide sense including unstable 
sphaleron solutions). 

Soliton solutions to non-Abelian Born-Infeld theory were discussed recently in a num- 
ber of papers using both the square root action with the ordinary and the symmetrized 



trace |[T0| , [LI], y, [12], [13[ , in the latter case, however, only perturbatively. It was shown that 
in the Born-Infeld theory, apart from monopoles and instantons, which also exist in the 
theory with the usual Yang- Mills quadratic action, new particle-like solutions of the glue- 
ball type are brought to existence |14| . For these solutions the full non-linear structure of 
the NBI lagrangian is essential, so earlier they could be explored only within the model 
with the ordinary trace. Here we show that these solutions persist in the theory with 
the action (|5]). Although the explicit Str lagrangian looks very differently from that with 
the ordinary tr, the glueball solutions remain qualiatively the same. We also study the 
static spherically symmetric NBI system coupled to gravity and demonstrate that gravity 
further reduces the difference between solutions obtained in these two models. 



2 Symmetrized trace NBI action for static 50(3)- 
symmetric fields 

We define the NBI Lagrangian in four spacetime dimensions as 



Lnbi = Str ^1 - ^ - det ( 9fiu + ^F)J) j = ^ Str(l - 11), 



(7) 



where 
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with the parameter (3 of the dimension of length 2 (the 'critical field'). The normalization 
of the gauge group generators is chosen as follows 



Ffiv — -^/iz/^aj tr t a tb — Sab- (9) 

The symmetrized trace of the product of p matrices is defined as 

Str(t ai . . . t ap ) = — Str (t ai . . . t ap + all permutations) , (10) 

and it is understood that the general matrix function like (|7|) has to be series expanded. 
It has to be noted that while under the Str operation the generators obviously can be 
treated as commuting objects, the gauge algebra can not be applied, i.e. r 2 7^ 1 until the 
symmetrization of the expansion is completed. 

A general 50(3) symmetric SU(2) gauge field is described by the Witten's ansatz 

V2A = a ti dt + aibi dr + {w t 2 - (1 - w) t 3 } dQ + {(1 - w) t 2 + w t 3 }sin6» d(/>, (11) 

where the functions ao, ai, w, w depend on r, t and v2 is introduced to maintain the 
standard normalization. Here we use a rotated basis t iy i = 1, 2, 3 for the SU{2) generators 
defined as 

t x = n a T a /y/2, t 2 = d e t 1 , smOt^d^h, (12) 

where n a = (sin 9 cos tp, sin 9 sin ip, cos 9), and r a are the Pauli matrices. These generators 
obey the commutation relations [U,tj] = -^-ijk^k- 

From four functions entering this ansatz one can be gauged away. In the static case 
we can further reduce the number of independent functions to two, while static purely 
magnetic configurations can be described by a single variable w(r). Here we deal with 
this simplest case 

y/2Ao = -(1 - w)t 3 , V2A V = sin 0(1 - w)t 2 , A r = A t = 0. (13) 

The field strength tensor has the following non-zero components 

V2F r9 = w%, V2F rip = - sin 9w't 2 , V2F 9v = sin 9{w 2 (14) 

where prime denotes derivatives with respect to r. For purely magnetic configurations 
the second term under the square root is zero, while the substitution of (|14D gives 

^9 , (1 — W 2 ) 2 9 W' 2 . 9 o. 

^ 2 = 1 + L ^V^^ + ^2 + ^)- (15) 

Now we have to expand the square root in a triple series in terms of the even powers of 
generators t\,t 2 ,tz- This can be achieved applying the formula 

(2m -2)! / x 



VTTx = \ ( TTT -7 > ( 16 ) 

^ 1 ml(m - 1)1 \ 4/ 



4 



and further repeatedly using binomial expansions. Finally we obtain 

j? ~ (_i)'+J+*(2i + 2j + 2fr - 2)! ^^^nW , 

L ™ = S ( JL (*+i + *-DK«l*! W IT] '<*•>•*>• (17 > 

where the sum over all positive k subject to a condition i + j ' + k > 1 is understood 
and 

V = 2fPr* ' K = W ( } 

and the final factor is the symmetrized trace of the product of even powers of Pauli 
matrices: 

f(i,j, AO = Str (rf r 2 2 Vf) . (19) 
To compute f(i,j, k) explicitly one can derive the following recurrent relation: 

(2i + 2j + 2k){2i + 2j + 2k — l)f(i, j, k) = 2i(2i - l)/(i - 1, j, k) + 

2j(2j-l)f(i,j-l,k) + 2k(2k — l)f(i,j, k — 1). (20) 

When only one index is non-zero one easily finds 

/(i,0,0) = /(0,i,0) = /(0,0,i) = 2. (21) 

The full solution of the Eq.([|), satisfying the border conditions ( pi]) reads 

_ 2(2Q!(2j)!(2ifc)!(» + j + ifc)! 
/( '' J ' fcj " (2< + 2j + 2ib)!iy!ib! ' ( j 



Substituting this explicit expression into the expansion (JL7\) of the lagrangian one obtains 
the following representation as a triple series 

/3 2 ^ (-l) i+ J +fc (2»)!(2j)!(2fc)! ^ 2 * /ifN^+ 2fc 
W " 4^ i!W(2i + 2j + 2k-l) [2) [JJ ■ {2S} 

Remarkably, one can perform this summation explicitly. First we observe that, once the 
factor l/(2i + 2j + 2k — 1) is omitted, and the summation is extended to all values of 
power indices including zero, we obtain a triple series expansion for the function 

Z(V, K) = 1 . (24) 

Indeed, treating this function as a product of three square roots, two of which are equal, 
we find the following representation: 



i,j,k=0 



Now it is easy to see that the desired sum ( p3|) is related to Z(V, K) through the following 
differential equation 

K^L(K, V) + V^L(K, V) - L(K, V) = £-{Z(V, K) — 1), (26) 



where L(K, V) stands for L NBI . Its solution satisfying the initial conditions 



L(0,0) = 0, _L(0,0) = — L(0) = 0, (27) 



following from an initial definition of the lagrangian (|7p, reads 



Lnbi = £ ( t _ l±£+£A , (2f 



where 



1 + V 2 V 2 - if 2 

^ = V^3^ arctanh VTTv^- (29) 

Here we assumed that V 2 > K 2 , otherwise an arctan form is more appropriate. Note that 
when the difference V 2 — K 2 changes sign, the function A remains real valued. It can be 
checked that when (3 — > oo, the standard Yang-Mills lagrangian (restricted to monopole 
ansatz) is recovered. In the strong field region our expression differs essentially from the 
square root /ordinary trace lagrangian. 



3 Glueballs 



The standard Yang-Mills theory does not admit classical particle-like solutions 1T3, 116 



17]. More precisely, this famous no-go result asserts that there exist no finite-energy 
nonsingular solutions to the four-dimensional Yang-Mills equations which would be either 
static, or non-radiating time-dependent fl7fl . This result follows from the conformal 
invariance of the Yang-Mills theory with the quadratic action, which implies that the 
stress-energy tensor is traceless : T£ = = — T 00 + Tu, where fi = 0,...,3, i = 1,2,3. 
Since T 0Q > 0, the sum of the principal pressures Tu is everywhere positive, i.e. the Yang- 
Mills matter is repulsive. This makes the mechanical equilibrium impossible JT8|. In the 



spontaneously broken gauge theories conformal invariance is also broken by scalar fields. 
Thus the above obstruction is removed what opens the possibility of particle-like solutions: 
magnetic monopoles and sphalerons. Monopoles are essentially related to the presence of 
the Higgs field with the SO (3) component, while for sphalerons Higgs can be replaced by 
other attractive agent. Recall that the sphaleron was first obtained in the gauge theory 
with doublet Higgs [|HJ and its existence was explained by Manton |21J as a consequence 



of non-triviality of the third homotopy group of the broken phase manifold. Later it was 
found that similar solutions arise in the theories without Higgs like Einstein- Yang-Mills 
and Yang-Mills with dilaton (for a review see [[H]] ) , in all such cases conformal invariance 
is broken. Recently it was observed that the Born-Infeld modification of the Yang-Mills 



action also breaks conformal invariance |TJ]], and gives rise to particle-like solutions. They 
form a discrete sequence labeled by the number of nodes of the function w(r), and the 
lower one-node solution is similar to the sphaleron of the Weinberg-Salam theory. 



In |14| the NBI lagrangian was adopted in the 'square-root '-ordinary trace form. Now 
we are able to perform similar calculations in the Str version of the NBI theory. It is 
worth noting that to study particle-like solutions in the pure NBI theory without Higgs 
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an exact in a' form of the lagrangian is needed since these solutions are formed in the 
strong field region. Therefore we assume the following one-dimensional action 



Note that the rescaling y/J3r — > r does not change the action (an overall factor appearing 
in the lagrangian can be removed by the corresponding rescaling of time). Therefore 
without loss of generality (3 can be fixed, it is convenient to choose (3 = 1/ y/2. Then the 
equation of motion for w will read 



w 



'K 2 VT+V* (W 2 -K 2 )A\\ Vw(K 2 A-V< 



2(V 2 -K 2 )\ 1 + K 2 ^/V 2 - K 2 )) (V 2 - K^VT+V 2 ' 



(31) 



where now 



y2 = ^M)! A" 2 — —pp- (32) 
Analyzing the extremal points of w as discussed in [ET| one finds that w can not have local 



minima for < w < 1, w < — 1 and can not have local maxima for — 1 < w < 0, w > 1. 
Thus any finite energy solution which starts at the origin on the interval — 1 < w < 1 lies 
entirely within the strip — 1 < w < 1. Once w leaves the strip, it diverges within a finite 
distance. We are interested in particle-like solution with finite total energy (mass) given 
for the present model by 



CO 



Boundary conditions at the origin can be derived combining the equation of motion with 
the requirement of convergence of this integral. Two classes of solutions are possible: one 
with w(0) = (leading to embedded U(l) solutions) and |w(0)| = 1, relevant for glueballs. 
Choosing without loss of generality w(0) = 1, we find the following series solution near 
the origin 

w = i _ br 2 + 36 2 r 4 326^+20^+3 6 

2(326 4 + 206 2 + 15) v ; v ; 

To ensure convergence of the total energy ( |33"D at infinity, w must tend to ±1, 0. In the 
case w(oo) = the global solution is w = 0, i.e. an embedded Abelian. For non-Abelian 
solutions one has 

w = ±(l-*)+o(\), (35) 



where c is another free parameter. The proof of existence of global solutions starting at 



the origin as (|34f) and approaching (|35|) at infinity may be given along the lines of |fL4 
The numerical integration shows that the discrete family of regular solutions exists for 
which the parameter b takes values shown in Tab. [1]. 

The integer n is equal to the number of zeroes of w. The n = 1 solution is an analog 
of the sphaleron known in the Weinberg-Salam theory |19|, 120], it is expected to have 



one decay mode. Higher odd-n solutions may be interpreted as excited sphalerons, they 
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n 


b 


M 


1 


1.23736 x 10 2 


1.20240 


2 


5.05665 x 10 3 


1.234583 


3 


1.67739 x 10 5 


1.235979 


4 


7.11885 x 10 6 


1.236046 


5 


4.94499 x 10 8 


1.2360497 


6 


4.52769 x 10 10 


1.2360497 



Table 1: Values of b and M for first six solutions 

are expected to have n decay directions in the configuration space. Even-n solutions 
are topologically trivial, they can be regarded as sphaleronic excitation of the vacuum. 
Qualitatively picture is the same as for the 'square root - ordinary trace form' but 



the quantized values of b are rather different in particular, in that case b% = 12.7463. 

Numerical solutions are shown on Fig. [I], where for comparison first four solutions 
with ordinary trace are also shown. The difference for the main n = 1 sphaleron in both 
theories is rather small, it increases for higher-ra solutions which move closer to the origin 
(where two models differ substantially). 

4 Gravitating glueballs and black holes 

Now let us consider the Str NBI theory coupled to gravity. The corresponding action 
reads: 

Senbi = J {RV=g + ^G(3 2 (Str J- det(^ + f3^F^) - y/=^ } d*x, (36) 

where R is the scalar curvature and G is the Newton constant. For static spherically 
symmetric solutions the metric can be parameterized as follows: 

dv^ 

ds 2 = Na 2 dt 2 - — - r 2 {d6 2 + sin 2 edxj>). (37) 



A computation of the symmetrized trace for the action ( |3"6D is a straightforward general- 
ization of the above procedure for the flat spacetime, the main difference being in using 
curved metric ( |3"TD instead of flat. After suitable rescaling, two dimensional parameters 
of the theory G, (3 combine in one dimensionless coupling constant 

9 = G/3, (38) 

which is the only substantial parameter of the theory. The reduced one dimensional action 
reads 
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where now 

and A is still denned by ( p9|) . 

The equations of motion derived from this action consist of an equation for the metric 
function a: 



&_ _ g K 2 r ((2V 2 + 2 K 2 V 2 - K 2 - K 4 )A- K 2 - K 2 V 2 ) 
~a ~ + V 2 {V 2 + K 2 V 2 -K 2 - K 4 ) ' 

an equation for the local mass function m(r) defined via N = 1 — 2m/ V: 

2 (l + V 2 + AK 2 



(41) 




m' = gr 2 y 7 === 1), (12) 

and the following equation for w: 



Nam' fK 2 Vl + V 2 (2V 2 - K 2 )A\ \ aVw(K 2 A-V 2 



2{V 2 -K 2 )\ 1 + K 2 Vl + V 2 Jj (V 2 - K^VT+V 2 ' 



(43) 



The equation for a decouples form the rest of the system, and this function can be 
found by a simple integration once N and w are known. Therefore we concentrate on a 
coupled system for N, w which is obtained after using the a-equation in the w-equation. 
Like in the ordinary trace version, these equations admit an embedded U(l) solution 
a = 1, w = corresponding to the unit magnetic charge. The corresponding metric was 
given in |22], [2| (see also @): 



N 



1 - \ ( m o + 9 J* (Vl + x 4 - x 2 ) dx^j , (44) 



where mo is a free parameter which can be positive, zero or negative. For mo = the 
black hole solutions (with the event horizon) exist for g > g cr = 1/2, otherwise there is no 
horizon. The role of this critical value in the non-Abelian case was discussed in [2~3| and 



(for the ordinary trace 577(2) Born-Infeld action). The metric for regular solutions 
approaches that of an Abelian solution (without horizon and mo = 0) for g < g cr and 
large n, the limiting mass being the corresponding Abelian mass. 

The present model also have regular gravitating solutions for all values of g which can 
be thought as interpolating between the flat case solutions discussed above and Bartnik- 
McKinnon |2£| solutions of Einstein-quadratic Yang-Mills model (for a detailed discussion 
see [0]). The situation is very much alike to the case of the ordinary trace ENBI model 



Regular gravitating solutions start at the origin with the following series expansion 



w = l-br 2 + 0(r^ 



3 VI + 46 2 V ; K J 
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where b is again a free parameter. At infinity one should have 

2M , s 

w^±l ,N~1 , 46 

r 

where M is the Schwarzschild (ADM) mass. Numerical solutions interpolating between 
these asymptotics are shown in Figs. |||. Note that the difference between Str and tr 
solutions is decreased with respect to the flat space case, especially for the first n = 1 
sphaleron. The metric bending is slightly more pronounced in the Str case. In Fig. f| we 
show the dependence of the parameter 61 for the n = 1 solution on the gravitational cou- 
pling constant for both Str and tr versions of the NBI theory. Qualitatively their behavior 
is the same. For g — > 00 both curves converge to the rescaled Bartnik-McKinnon's value. 

In weak gravity limit (as g — > 0), u>(r)-functions for regular solutions do not differ 
considerably from the flat case, especially the sequence b n is unbounded while the number 
n of nodes of w tends to infinity. The metric function N(t) with increasing n approaches 
the metric of the abelian solution at some interval which moves more and more close to 
the origin. The ADM masses of these solutions behave like 

]i m MM=Ml lat (47) 
9-*o g 

where M^ at are the flat case masses defined by (|33| ) and shown in table 1. 

But if g > gw = 1/2 the situation changes. With increasing n the parameters b n tend 
to a limiting value and the metric functions tend to the metric of the limiting abelian 
solution with a degenerate horizon. This situation resembles that of the Einstein- Yang- 
Mills model and indeed this model could be recovered in the limit g — > 00 after a rescaling 
^fgr — > r. In this limit the glueball masses behave like 

lim MM = m bk (48) 

where M^ K is the mass of corresponding Bartnick-McKinnon solution for the Einstein- 
Yang-Mills model. This behavior could be observed in Fig [5]. 

Let us discuss the black holes. Instead of ( |45|) now we specify boundary conditions at 
the horizon r = using the following series expansions: 

N = N' h (r-r h )+0((r-r h ) 2 ), 

*> = w h+^kA(r-r h ) + 0((r-r h ) 2 ), (49) 

where 



N' h = -2gr 2 h (VTTK^ - 1) + 1. (50) 
r h 

Asymptotically flat black hole solutions exist for all horizon radii rh and some discrete 
sequence of Wh, labeled by the number of nodes of w like in the regular case. Qualitative 
behavior of w outside the horizon is also the same. Numerical results for the metric 
function N(r) of the n — 1,3 black holes are shown in Fig. ^ for a critical g = 1/2 and 
different values of the horizon radius. For non-small the function N is monotonous 
outside the horizon, while for smaller one observes a local minimum near the black 
hole. 
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5 Discussion 



Two main results of this paper should be noted. First, we obtained a closed analyti- 
cal expression for the symmetrized trace version of the SU(2) NBI action restricted to 
monopole ansatz. This expression is particularly simple and can be further used in various 
problems. Moreover, similar technique can be applied to other related field configurations, 
the results will be given elsewhere. 

Second, we have extended results of jnj about the existence of sphaleronic glueballs 
in the 577(2) NBI theory without Higgs field to the symmetrized trace version of the NBI 
model. In both cases the classical scale invariance of the ordinary Yang-Mills lagrangian 
is broken by the Born-Infeld non-linearity, what removes an obstruction for classical glue- 
balls. We have found that qualitatively the Str NBI glueball solutions remain the same, 
although a certain difference is observed in the core. In particular, quantized values of 
the parameter determining the derivative of the Yang-Mills field at the origin are much 
larger in the Str case. Nevertheless, the masses are not very different and with growing 
node number converge rapidly to the same limit. This is due to the fact that most of the 
energy is localized in the outer core region where the difference between two models is 
less pronounced. 

When gravity is taken into account, we have demonstrated that there exist a con- 
tinuous transition between the flat space glueballs and Bartnik-McKinnon's particle-like 
solutions of the Einstein- Yang-Mills theory with the usual quadratic lagrangian. The lat- 
ter can therefore be viewed as the strong gravity limit of purely gluonic sphalerons on 
D3 branes. Gravity was shown to further reduce the difference between the results of the 
Str and tr models. This is well understood, since for the strong gravity (in units of the 
Born-Infeld critical field) both models have the same limit. 

We thank A. Koshelev for a useful discussion. This work was supported in part by 
the RFBR grant 00-02-16306. 
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Figure 1: First six solutions w n (r) for flat space glueballs in the NBI theory with sym- 
metrized trace (solid line) and their four ordinary trace counterparts (dotted line). 
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Figure 2: Functions w n (r), n — 1, 2, 3 for gravitating glueballs in the Str (solid line) and 
Tr (dotted line) NBI models with g = g cr = |. Solutions are practically indistinguishable 
for n — 1, becoming slightly different for n = 2,3. 
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Figure 3: Metric functions N n (r), n = 1,2,3 for gravitating glueballs in the Str (solid 
line) and Tr (dotted line) NBI models with g = 1/2. Gravitational binding in the Str case 
is slightly greater than in the Tr case. 
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Figure 4: Parameter 6 versus an effective coupling g (in log variables) for the gravitating 
NBI model n — 1 solutions in the Str (solid line) and Tr (dotted line) models. 
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Figure 5: Dependence of mass on g for the n = 1 gravitating glueball (in units y/g) in the 
ENBI theory with symmetrized trace (solid line) and ordinary trace (dotted line). 
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